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1 . I nt roduct ion 

Most investigations of the muon arrival time distribution in EAS during 
the past decade made use of parameters which can collectively be called 
rise time parameters. We will follow the definition of Blake et al. 

(1981) and define the rise time parameter T^/g as the time taken for the 
integrated pulse from a detector to rise from Al to Bl of its full 
amplitude. The use of these parameters are usually restricted to the 
determination of the radial dependence thereof. This radial dependence 
of the rise time parameters are usually taken as a signature of the par- 
ticle interaction characteristics in the shower. As these parameters 
have a stochastic nature, it seems reasonable to us that one should also 
take notice of this aspect of the rise time parameters. The aim of this 
paper is therefore to present a statistical approach to the rise time 
parameters, as this has not been done in the past. 

2 . Order statistics and rise time parameters 

From the definition of T^/g it is reasonable to assume that on the ave- 
rage Al of the total number of particles which gave rise to the output 
pulse from the detector have arrived in the time Tq/a. A sample quantile 
of order p is defined as follows: Let (X-j, X2 . ..X n ) be a random sample 

of size n of a random variable X with probability density function (pdf) 
f(t). Let (X| n | , X£ n ' . . .X^h) ) be an arrangement of (X^, X2 ... X n ) such 

that xj"> £ xl n £ •*.X^ n ^. x| n ) iscalled the i-th order statistic. The 
sample quantile of order p, 0<p<1, is defined to be the order statistic 
X|< f° r which k = [np] + 1, where [np] is the greatest integer not larger 
than np. This definition states that a fraction p of the sample values 
is less than x£ n ' . This corresponds exactly to the definition of Tq^. 
The rise time parameter T^g corresponds therefore to the difference 
between two order statistics of the sample. The statistical properties 
of T^/g can therefore be determined from the properties of the difference 
between two order statistics. 


We now give a number, of properties which are very useful. The following 
notation is used: fA n ' (t) is the pdf of X£ n ' , Wj^(t) is the pdf of 

X} n '-X£ n ', F(t) is tne distribution function of the random variable 
X, i.e. F ( t ) = P(X<t), and g n '(t) is the pdf of the arrival times of 
muons with respect to the first detected muon. 


(i) 

(N) 


f ( n ) ( t ) .r’ C f(t)[F(t)] k 1 [l -F(t)] n k with C = n ! / (k- 1 ) ! (n-k) ! 

^ 00 

W, {n J(t) = C [ f (t)f (x+t) [F(t)] k 1 [ F (x+t) -F (x )] ^ k 1 [l-F(x+t)] n k dx 


'o 

with C = n !/(k-l) ! (n-k) ! (€-k-1 ) ! 
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(Swanepoel (1985)). 


Property (Ml) states that the sample quantile converges stochastically to 
the population quantile. Properties (iy). (v) and (vi) give the limiting 
distributions of f£ n '(t), Wjjjp (t) and g(")(t). The rise time parameters 
will then also have the above properties. 

The significance of these results is further that it is possible to examine 
specific properties of the rise time parameters without having to simulate 
air showers. It must be made clear at this point that points (i) to (vi) 
does not add directly to the understanding of the physics involved in 
shower development. We feel, however, that it gives one a firm basis from 
which one can understand the properties of the measured rise time para- 
meters. 

Discussion 


Consider, for example, property (iii). For the case of EAS, the arrival of 
the hypotherical shower front can be considered as the time t=0 from which 
timing measurements can be made. In this case f(t) = 0 for all t < 0, and 
we can allow p to be equal to zero. Then Xj n J, corresponds to the~zeroth 
order sample quantile while a = a D = 0. We then have 1 im P( |x| n '-0 |>e) 

= 0, which means that the arrival time of the first detected muon with 
respect to the shower front converges stochastically to zero. This means 
that for large samples one can consider the first detected particle as the 
time t - 0. Property (vi) is equivalent to what has just been said. 

In figure 1 we show the arrival time distribution of muons together with 
the arrival time distribution of muons when the first detected muon was 
taken as the time t = 0. T he difference at small delays is due to the fact 
that the sample size was only ten. 

In figure 2 we present the relationship between <X>5)- x|5)> anc | <y> for the 
case when f(t) is a gamma density function. The same relationship was 
found to exist for larger samples. Is is also possible to show that this 
linear relationship exists for a wide variety of density functions. In 
figure 3 we present the relationship between <xjjj '- X'j°'> and <X> for 
shower simulation data. Xp®) is the i-th order statistic for a sample 
of 10 muon arrival times measured in a detector and <X> is the mean muon 
arrival time. It can be seen that the relationship is also linear even 
though we do not know the parametric form of the muon arrival time distri- 
bution. This example also illustrates that it is not necessary to 
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determine the properties of the rise time parameters through the simula- 
tion. of air showers. It should also be noted from figure 2 that 
<X^ n '- X|| n '>does not determine <X> uniquely. This is an example of 
some of the limitations of rise time parameters. 

Conclusion 


With the above examples we have tried to illustrate that there exist a 
statistical basis for the rise time parameters. We bel ievp that the 
statistical properties of rise time parameters may be of use not only for 
the analysis of experimental data but also for the planning of experiments. 
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Figure 1 Comparison between g ^(t) and f(t) 

for 14 Fe- initiated showers. f(t), 

— _ g ( 1 0) ( t ) < 110m < R < 120m, 

E y > 1 GeV 
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